Adiabatic Quantum Transport in a Spin Chain with a Moving Potential 
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Many schemes to realize quantum state transfer in spin chains are not robust to random fluc- 
tuations in the spin-spin coupling strength. In efforts to achieve robust quantum state transfer, 
an adiabatic quantum population transfer scheme is proposed in this study. The proposed scheme 
makes use of a slowly moving external parabolic potential and is qualitatively explained in terms 
of the adiabatic following of a quantum state with a moving separatrix structure in the classical 
phase space of a pendulum analogy. Detailed aspects of our adiabatic population transfer scheme, 
including its robustness, is studied computationally. Applications of our adiabatic scheme in quan- 
tum information transfer are also discussed, with emphasis placed on the usage of a dual spin chain 
to encode quantum phases. The results should also be useful for the control of electron tunneling 
in an array of quantum dots. 
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I. INTRODUCTION 

Efficiently transferring quantum information is an im- 
portant challenge for the practical realization of quantum 
computers. Optical fibers, where quantum information 
is transmitted by mobile carriers like photons, are most 
desirable for long-distance communication. However, the 
need to interface with solid state quantum computer com- 
ponents considerably restricts the experimental feasibil- 
ity of using optical systems as quantum channels. One 
promising alternative is to use condensed matter systems 
themselves as the quantum wire for information transfer. 
This idea, initially advocated and studied by Bose [ij in 
the context of quantum spin chains, has now attracted 
wide interests. In particular, for a quantum spin chain 
used as the quantum wire, the natural evolution of per- 
manently coupled spins is exploited to accomplish the 
quantum information transfer. Because the Hamiltoni- 
ans equivalent to that of a spin chain may be realized 
in different physical systems (e.g., arrays of Josephson 
junctions cold atoms in optical lattices arrays of 
quantum dots, etc.), interests in spin chains as a promis- 
ing candidate of quantum wires continue to grow. Exper- 
imental studies on the dynamics of quantum information 
transfer and on the quantum control in spin chains us- 
ing nuclear magnetic resonance techniques have also been 
reported Q. 

Despite many fruitful studies of spin chains from 
a quantum information perspective, many theoretical 
problems are still open. To be more specific let us con- 
sider first Bose's original proposal [l| . Therein the fidelity 
of quantum information transfer is gradually degraded by 
the dispersion effects associated with the quantum prop- 
agation. Furthermore, particular external magnetic field 
should be designed to ensure a correct quantum phase at 
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the receiver's site. These issues and others motivated a 
series of sophisticated protocols to achieve better quan- 
tum information transfer. One noteworthy approach was 
to pre-engineer the nearest-neighbor couplings of a spin 
chain or even a spin network 0, || @, H, A sec- 
ond approach exploits the mirror symmetry of a spin 
chain [H, [l(| ■ Another approach s ugg ests to use Gaussian 
wavepackets with slow dispersion [ll| to encode the quan- 
tum information to be transferred along an unmodulated 
spin chain. Unfortunately, these pioneering approaches 
rely upon specific analytical forms of the involved Hamil- 
tonian and hence are not robust to imperfection or phys- 
ical fluctuations in the spin-spin coupling strength. One 
encounters the same situation when applying other more 
subtle techniques [llGjJEi]. 

To have the desired robustness that may be necessary 
for any type of quantum information transfer, two quan- 
tum control schemes based on adiabatically varying the 
cou plin g strength in a spin chain have also been suggested 
fl5lll6|]. Note however, these schemes require individual 
addressing of the nearest-neighbor coupling and hence 
present new experimental challenges. Another novel and 
quite robust quantum transfer protocol in spin chains is 
the dual spin chain scheme [13, EH Ell H3| • Therein the 
quantum information to be transferred is encoded in two 
parallel spin sub-chains (initially assumed to be identi- 
cal, but even this condition may be lifted under certain 
conditions). Thanks to the encoding with two spin sub- 
chains, the quantum transfer can be very robust to static 
disorder. Nevertheless, even this promising scheme is not 
perfect, because (i) it may need too many quantum mea- 
surements (or too many steps of "trial and error"), (ii) it 
may not operate well in the presence of time-dependent 
disorder, and (iii) the effects of nonideal measurements 
are still under investigation fl9j | . 

Hence, it remains an open question as to which quan- 
tum transfer scheme will ultimately be adopted experi- 
mentally, with high fidelity and low cost. It is our belief 
that in the end a combination of several techniques may 
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be able to offer the most powerful protocol for quantum 
information transfer in solid state systems. 

In this paper, we introduce an adiabatic transport 
scheme assisted by a slowly moving external field applied 
to a spin chain. Thanks to a pendulum analogy, the cen- 
tral idea can be understood in a very simple manner. The 
essence is that when an external field is moving slowly, 
the spin excitation may adiabatically follow the field un- 
der certain conditions. During this process the quantum 
population of spin excitation is transferred from one end 
of the spin chain to the other end. As we show below 
in detail, this adiabatic scheme offers a number of ad- 
vantages: (i) it is highly robust, (ii) it can be operated 
rather fast (in the absence of disorder, it may be as fast 
as one tenth of the natural propagation speed of the spin 
chain), (iii) the required field strength can be decreased 
by using wavepackets as initial states, (iv) the transfer 
can be easily stopped and relaunched, (v) the time of ar- 
rival of the quantum population transfer to the last spin 
with high probability can be easily predicted, and (vi) 
it can offer a promising means of quantum information 
transfer when combined with the above-mentioned dual 
spin chain scheme. 

The outline of the paper is as follows. In Sec. [Il] we 
explain the motivation and mechanism of our adiabatic 
population transfer scheme by mapping the spin chain 
Hamiltonian to that of a pendulum. Detailed computa- 
tional results are presented in Sec. III. The robustness 
of our scheme is studied in Sec. IV. In Sec. V, we dis- 
cuss how our adiabatic population transfer scheme, which 
does not yet take care of the quantum phases (also essen- 
tial for quantum information transfer), can be combined 
with the dual spin chain scheme to offer a potentially 
powerful approach for quantum information transfer. We 
conclude this work in Sec. VI. 



II. ADIABATIC QUANTUM TRANSPORT IN 
SPIN CHAINS: A PENDULUM PERSPECTIVE 

Consider a one-dimensional Heisenberg chain of N + 1 
spins subject to an external parabolic magnetic field. The 
associated Hamiltonian is given by 

jN-x n c 

H s = ~-^a n -(T n+1 + ^—{n-n Q ) 2 a z n , (1) 

n=0 ri=0 

where a = (a x , a v ,a z ) are the Pauli matrices, J the cou- 
pling strength between nearest neighbor spins, and C is 
proportional to the magnetic dipole of the spins and the 
amplitude of the parabolic field whose minimum is at site 
no. Note that uq will be time-dependent in our control 
scheme. Below we assume all the system parameters have 
been appropriately scaled and take dimensionless values, 
with J = 1, h = 1 throughout. As such, the energy scale 
(e.g., the parameter C) should be understood with re- 
spect to J, and the time scale should be understood with 
respect to Ti/J. Because the spin chain Hamiltonian H s 



commutes with the total polarization S z = Xm=i a ni ^ ne 
dynamics of the spin chain preserves the total polariza- 
tion. Here we restrict our analysis to the subspace of 
S z = 1 — N, where in total only one spin is flipped. In 
this subspace the total state of the chain can be written 
as 

N 

!*(*)} = E c ™Wl m >> ( 2 ) 

m=0 

where |m) represents one basis state with a spin up at 
the mth site and all other spins down. The complex 
coefficients c m (t) are the probability amplitude. Below 
we also shift the zero of the energy scale such that if one 
spin is down, its interaction with the external magnetic 
field contributes zero to the total energy. 

To understand the essence of the spin chain dynam- 
ics from a semiclassical perspective, we now consider the 
large N limit of the spin chain. Denote k as the quasi- 
momentum of a plane spin wave and the \k) the associ- 
ated eigenstate of the quasi-momentum. Then the Hamil- 
tonian in Eq. {T]) can be rewritten as 

H s = -J cos{k)\k)(k\ 
Jo 

+ 5>)(n||(n-n ) 2 . (3) 

n 

This form can be further simplified in an operator form, 
i.e., 

C 

H s = -Jcos(k) + -{h-n ) 2 , (4) 

where k\k) — k\k); n\n) — n|n); [cos(fc),n] = — isin(fc), 
and [sin(fc),n] = icos(k). The Hamiltonian in Eq. (H} 
can now be easily recognized to be the Hamiltonian of 
a quantum pendulum H p (2lj . with an effective Planck 
constant \fC . Specifically, with the mapping k — > x, 
\J~Cn — > p, the spin chain Hamiltonian H s is mapped to 

H P = -Jcos(x) + i(p-po) 2 , (5) 

where po = vCno and [cos(x),p] = —i\/C sin(x). The 
semiclassical Hamiltonian for this quantum pendulum, 
i.e., Hp = — Jcos(x) + \{jp — po) 2 is obtained by replacing 
x and p with c- variables x and p. 

Many quantum transport features of the spin chain can 
now be understood in terms of the semiclassical dynamics 
of the pendulum analogy thus obtained. In particular, 
the quantum transport in the momentum space of the 
pendulum is now in parallel with the transfer of spin 
excitation from one site to another. Hence, the issue 
of robust quantum transport of spin excitation along the 
spin chain now reduces to the design of a control scenario 
that enables robust transport of the pendulum state in 
its momentum space. 
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Schrodinger equation, 




FIG. 1: Phase space portrait for the classical pendulum 
Hamiltonian associated with Eq. ((5j), with C = 2 and J = 1. 
The closed curve enclosing the island and intersecting with 
both p — po — and x = ±7r is the separatrix. Variables 
P — Po and x plotted here take dimensionless values. Note 
that if the separatrix is moving up slowly by increasing po 
gradually, an initial state enclosed by the separatrix is ex- 
pected to adiabatically follow the movement of the separa- 
trix. This suggests that a slowly moving parabolic potential 
applied to a spin chain can be used to adiabatically transfer 
spin excitation along a spin chain. 



An adiabatic scheme for robust quantum population 
transfer along a spin chain can now be proposed. The key 
observation is the existence of a motional separatrix in 
the classical phase space of the pendulum. This separa- 
trix is located at — Jcos(a;) -I- (1/2) (p — p ) 2 = J [see Fig. 
[Tj. If we now slowly move up the separatrix along the 
momentum space by increasing po, then a quantum state 
initially trapped inside the separatrix cannot penetrate 
this separatrix and is expected to adiabatically follow the 
moving separatrix, giving rise to adiabatic transport in 
the momentum space. Translating this pendulum lan- 
guage back to the spin chain case, one anticipates that a 
slowly moving parabolic magnetic field (with slowly in- 
creasing no) should result in a robust scenario for trans- 
ferring quantum population along the spin chain. During 
this process the dispersion of the spin wave should also be 
bounded by the separatrix structure, i.e., a moving but 
non-spreading wavepacket [22| of the spin wave can be ex- 
pected. The main remaining task of this paper is devoted 
to detailed aspects of this adiabatic control scheme. 



It is also interesting to note that the dynamics of 
the spin chain can be mapped to that of a tight- 
binding model. To see this consider first the associated 
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(c„_i + Cn+i) + —(n - n ) On, < n < N, 
C 

Cjv-i + TT {N - n Q ) 2 c N . (6) 



Consider next a tight-binding Hamiltonian H t describing, 
for example, an array of (N + 1) identical quantum dots 
subject to an external parabolic field, with one electron 
tunneling between the quantum dots. Then H t assumes 
the following form, 



H t = --x J2(aia n+1 +a n a ] n+l ) 



N-l 



N 



c 



— (n - n ) 2 4a„, 



(7) 



n=0 



where J represents the constant tunneling rate between 
the nearest-neighbor quantum dots, and and a n rep- 
resent the creation and annihilation operators. Because 
the total number of electrons is already assumed to be 
one, the system wavefunction can also be written as 

= Em=o c m(0l m )> where |m) denotes the state 
with an electron in the mth quantum dot and c m (t) de- 
notes the associated quantum amplitude. In this rep- 
resentation, one immediately finds that the evolution of 
this tight-binding system takes the same form as Eq. ^ . 
Thus, the above pendulum analogy is also applicable to a 
tight-binding system and is hence very useful for consid- 
eration of adiabatic quantum transport in quantum dot 
arrays 23] . Note also that one may start from Eq. ((6|) to 
have an alternative derivation of the pendulum analogy 

To end this section, we stress that the proposed control 
scheme is based upon a semiclassical perspective afforded 
by the pendulum analogy. What is not addressed is the 
issue of transferring the quantum phase along the spin 
chain. As such, although the pendulum analogy helps 
design our scheme for the robust transport of quantum 
excitation along a spin chain or the robust transport of an 
electron in a quantum dot array, the issue of quantum in- 
formation transfer is only partially touched. Indeed, the 
introduction of an external field will change the energy 
of the spin chain system and hence will necessarily in- 
troduce extra dynamical phases to the evolving quantum 
system. This makes it clear that transporting quantum 
phases along the spin chain requires additional consider- 
ations. This quantum phase issue will be considered in 
detail in Sec. V. 



III. ADIABATIC TRANSPORT BY A MOVING 
POTENTIAL: COMPUTATIONAL RESULTS 

In this section, we illustrate our adiabatic quantum 
population transfer scheme with detailed computational 
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FIG. 2: Excitation probability transferred to the last spin in 
adiabatic quantum transport along a chain of 101 spins, as 
a function of the amplitude of the external moving magnetic 
potential characterized by C [in units of J, see the text below 
Eq. (1)]. The moving speed of the control field is S = 0.005. 

examples. Let us assume that the initial state of the spin 
chain is given by |$) = J2m=o c »n(0)l m )- Two types of 
Cjn(0) will be considered below. In the first case only 
the m = 0th spin is excited, with c m (0) = 5 m Q. In the 
second case, the initial state is a Gaussian wavepacket 
truncated to three sites only, with c m (0) cx exp[— (m — 
l) 2 /2/o] for m = - 2 and l = 0.707 being the width 
of the Gaussian wavepacket. In either case a parabolic 
magnetic field first centered on the n = 0th site is applied 
and then slowly moved to the regime of larger n. This 
is realized by introducing the time dependence of no via 
n = + St, where S is the moving speed. Note that a 
static parabolic field was previously introduced to induce 
a quasi harmonic lower energy spectrum such that good 
transfer of Gaussian wavepackets [l0( may be realized. 
By contrast, our moving potential scenario is more active 
and effective in controlling the quantum transport and is 
in principle applicable to cases where the shape of the 
external potential is not parabolic. 

The quantum state of the spin chain at a later time 
can be directly calculated using the Schrodinger equation 
given above. In particular, the probability of transferring 
the quantum excitation to the last spin of the chain can 
be examined. If the performance of the population trans- 
fer is satisfactory, one should find \cm\ 2 ~ 1-0. Evidently, 
this condition of high transfer probability is already use- 
ful by itself for, e.g., transporting electrons in a quantum 
dot array in a controlled fashion. As shown in Sec. V, the 
phase of the quantum amplitude cjv may be also taken 
care of by considering a dual spin chain. 

We now discuss the feasibility of adiabatic quantum 
population transfer by taking advantage of the separatrix 
associated with the pendulum analogy. In the ideal case 
of adiabatic following, an initial quantum state enclosed 
by a separatrix will move with the slowly moving separa- 
trix. Consider first an initial state localized exclusively 



at the n — 0th site. Then the associated fc-distribution 
covers uniformly from to 2tt. From a semiclassical per- 
spective afforded by the pendulum analogy, such an ini- 
tial state corresponds to an initial ensemble lying on the 
(p — po) = axis of the classical phase space. This initial 
ensemble hence necessarily intersects with the separatrix 
(see Fig. [1]). Because the motional period associated 
with the separatrix is infinity, those ensemble compo- 
nents that overlap with the separatrix will always regard 
the movement of the separatrix as "too fast to follow" . 
That is, as we slowly move the separatrix upwards in the 
classical phase space, some portion of the initial ensemble 
may break the adiabaticity and tunnel through the sepa- 
ratrix structure. Under such a situation adiabatic quan- 
tum population transfer is expected to partially break 
down. To reduce the degree of non-adiabaticity, one pos- 
sible approach is to reduce the overlap of the initial state 
with the classical separatrix. This should be doable by in- 
creasing the effective Planck constant y/C (i.e., increasing 
the strength of the parabolic field) such that the separa- 
trix regime supports less quantum states. This is indeed 
what we find computationally for a chain of 101 spins. 
In particular, Fig. [2] shows that for a field amplitude 
characterized by C — 0.5, the probability of transferring 
the initial excitation to the last spin is only 0.63. By 
increasing C to 8.0, a transfer probability around 99% 
is observed. Figure 3 shows the actual excitation pro- 
file |c„| 2 vs. n for a spin chain subject to a parabolic 
potential moving at a constant speed of S — 0.005. At 
t = 0, the state is at site n = 0. At t = 10000, the 
quantum population is mainly at n = 50. Note that at 
that moment the excitation profile is slightly delocalized 
into three sites, but the peak probability is still as high 
as 0.97. This peak is propagated to site n — 100 at time 
t = 20000, with no further dispersion detected. This indi- 
cates that our moving potential scheme has the capacity 
to overcome the dispersion issue in quantum information 
transfer. Though the required field strength for large 
| n — no | could be demanding experimentally, we point 
out that because the spin excitation is highly localized 
throughout the process, the moving parabolic magnetic 
field does not need to span over many spin sites (in our 
numerical experiments, we use a parabolic field that only 
spans 20 sites). 

We have also examined the quantum dynamics for the 
second type of initial states, i.e., states with a Gaussian 
excitation profile at t = 0. Because such initial ensem- 
bles are localized in both k and n, they can be naturally 
enclosed by the separatrix shown in Fig. [TJ As such, if 
the shape of the initial excitation profile is appropriately 
adjusted, the initial state can be made not to intersect 
with the separatrix. This being the case, the adiabatic 
following should work better, probably requiring a much 
weaker parabolic field. This expectation is also confirmed 
computationally. In particular, Fig. 3] shows the trans- 
port of an initial Gaussian excitation profile, again for 
a chain of 101 spins. At t = 0, the excitation profile 
spans only the first three sites with a probability peak 



■5 



1.0r- 
0.8 
-9 0-6 

i " 

0.2 
0.0 



1.0r 



(a) 



o.s (b) 



■O 0.6 ■ 
S 0.4- 
0.2 

o.o L 



1.0 r 



0.8 (C) 
jd 0.6 - 
2 0.4- 
°- 0.2- 

0.0 1 



40 50 i 

Site index 



40 50 
Site index 



40 50 I 

Site index 



(a) 



10 20 30 40 50 60 70 
Site index 





10 20 30 40 50 60 70 80 90 100 
Site index 



FIG. 3: Adiabatic transfer of spin excitation initially local- 
ized exclusively at the n = 0th site along a chain of 101 spins, 
for a field amplitude given by C = 8 and its moving speed 
given by S — 0.005. Panels (a), (b), and (c) are for times 
t = 0, 10000, and 20000. Note that the final peak excitation 
probability remains as high as 0.97. 



0.78 at the site n = 1. This state is then transported by 
applying a parabolic field with C = 2 moving at a rate 
of S = 0.005. During the quantum transport the state 
disperses among about five sites, with a peak probability 
maintained around 0.77. At t = 10000 and t = 20000, 
the peak of the spin excitation probability profile is trans- 
ferred to n — 50 and n = 100. Interestingly, it is found 
that the population transfer probability to the last spin 
can be further enhanced at a slightly later time. As seen 
from Fig. [4j at time t = 21000, the peak probability, 
located at the last spin, is as high as 0.997. Physically, 
this is due to the reflection process at the end of the spin 
chain. In some sense, the interplay of the parabolic field 
centered at the end of the spin chain and the reflection 
process acts as a lens refocusing the slightly dispersed 
profile, and the peak probability builds up on the last 
spin. Note also that in the absence of the control field, 
one in general needs initial wavepackets of much larger 
length to be able to reduce the undesired dispersion [Tl| • 

Results in Fig. 2] demonstrate that using initial states 
whose excitation profile covers a few spins can signifi- 
cantly reduce the required field amplitude (compare val- 
ues of C in Fig. |4] and in Fig. [3]) . Counter-intuitively, 
there exists a maximal number n max of spins that can be 
used to create such initial states. This can be appreci- 
ated by considering again the separatrix in the classical 
phase space of the pendulum analogy. Classical orbits 
outside the separatrix are associated with pendulum's 
rotational motion (rather than oscillation). Going back 



FIG. 4: Adiabatic quantum transport along a chain of 101 
sites with an initial Gaussian excitation profile, at (a) t = 
0, (b) t = 10000, (c) t = 20000, and (d) t = 21000. The 
amplitude of the parabolic field is given by C = 2, and the 
moving speed of the control field is given by S = 0.005. Note 
that the final excitation probability transferred to the last 
spin is as high as 0.997. 



to the spin chain or the tight-binding system, these states 
correspond to Bloch oscillations in a "locally linear" field. 
As the separatrix is slowly moving, these states can con- 
tinue their Bloch oscillations in a slowly-varying local 
field and hence will not follow the motion of the separa- 
trix in the momentum space. With this understanding, 
one may estimate n max from the width of the separatrix 
in the momentum space. Specifically, for a fixed value of 

the parameter C, n max ~ + 1. For the numerical 

example in Fig. 2J one obtains n m ax ~ 4. This estimate is 
quite consistent with the finding that during the popula- 
tion transfer, the moving wavepacket does not cover more 
than five sites. This result also implies that for weaker 
magnetic fields (smaller C), one can use more spins to 
form the wavepacket for analogous adiabatic population 
transfer. 

Because our quantum transport scheme is based upon 
the adiabatic following of the spin excitation profile with 
a moving external potential, it can stop and relaunch 
the excitation transfer at any time with great ease, by 
simply stopping and restarting the movement of the ex- 
ternal parabolic potential. This is simpler than a recent 
approach [25j using pulsed magnetic fields, and is also 
confirmed in our computational studies (not shown). 

So can we further increase the moving speed of the con- 
trol potential while still maintaining the adiabatic follow- 
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FIG. 5: Adiabatic transfer of spin excitation for an initial 
state exclusively localized at the n = 0th site along a chain of 
101 spins. The amplitude of the moving parabolic potential is 
given by C = 8, and the moving speed is given by S = 0.025. 
Panels (a), (b) and (c) show the excitation profile at times 
t = 1000, 2600, and 4000. 



ing and hence the adiabatic quantum population trans- 
fer? Our findings in this regard can be summarized as 
follows: (i) for large C adiabatic quantum transport may 
survive for a moving speed around 10% of the coupling 
constant J. The smaller the field strength C is, the lower 
the threshold moving speed will be; (ii) when the moving 
speed exceeds the threshold, the probability of successful 
population transfer gradually decreases, but can still be 
considerably large for a relatively short spin chain. For 
example, Fig. [5] shows the result for an initial state ex- 
clusively localized at the n = 0th site. The moving speed 
of the parabolic potential is S = 0.025. The peak value of 
the probability profile is 0.96 at time t = 1000. It reduces 
to about 0.95 at t = 2600 and 0.94 at t = 4000. Hence, 
in this case, only 2% reduction in the peak probability 
occurs when the moving speed S increases by a factor 
of five. However, increasing the moving speed beyond 
this limit drastically reduces the probability of popula- 
tion transfer to the last spin. For a moving speed of 
S = 0.30, the probability maxima equals only 0.87 at 
t = 2200, and 0.84 at t = 3400. Analogous calculations 
are also carried out for the transport of an initial Gaus- 
sian excitation profile. As shown in Fig. [6] for a moving 
speed of 0.1, the peak value of probability remains around 
0.76 during the transport process. As such, at the end 
(not shown) the adiabatic population transfer is also very 
successful for this high moving speed. But if the moving 
speed is further increased by several times, dispersion in 
the spin excitation profile will be considerable. 

In short, our numerical experiments suggest that, to 



FIG. 6: Adiabatic transport of an initial Gaussian profile of 
spin excitation (same as in Fig. [4} along a chain of 101 spins. 
The amplitude of the moving parabolic potential is given by 
C = 2, and the moving speed is given by S — 0.1. Panels (a), 
(b), and (c) are for t = 300, 600 and 1000. 

achieve adiabatic transport of spin excitation along a 
quite long spin chain using a moving parabolic potential, 
the associated moving speed can be as large as being one 
tenth of the natural propagation rate (J) of the system 
(without disorder). 

IV. ROBUSTNESS OF ADIABATIC 
TRANSPORT 

So far the majority of quantum state transfer schemes 
consider only idealized spin chains with no disorder in the 
spin-spin coupling strength. This suggests a gap between 
theoretical exploration and realistic situations in exper- 
iments. In particular, the effects of static and dynamic 
imperfections in spin chains are studied in very few cases 
[13l Il8l . |20L l26j . Here we computationally study the influ- 
ence of static and dynamic disorder on our adiabatic pop- 
ulation transfer scheme, by considering the model Hamil- 
tonian in Eq. {1} with fluctuating spin-spin coupling con- 
stants. We hope to numerically confirm the robustness 
of our scheme as implied by its adiabatic nature. 

The model Hamiltonian with disorder in the spin-spin 
coupling strength is given by 

„ ^ (J + S n ) 

Hsd — 2 ' 11+1 

n=0 

N „ 

+E 2>- n °) 2 <' (8) 

n=0 

where S n are time-independent random numbers uni- 
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FIG. 7: Adiabatic transfer of spin excitation for an initial 
state exclusively localized at the n = 0th site along a chain 
of 101 spins, in the presence of static disorder with the noise 
amplitude given by A = 0.5. Other parameters are the same 
as in Fig. |3] The excitation profile at time t = 7000, 14000, 
and 20000 are shown in panels (a), (b), and (c). Results here 
are very similar to those shown in Fig. [3] 



formly distributed in the interval [—A, A], representing 
random fluctuations in J with the amplitude A. We call 
cases with such time-independent disorder as static dis- 
order models. Note that specific results presented below 
refer to single disorder realizations. These results are 
very typical such that there is no need to average over 
many disorder realizations. 

Figure [7] display one sampling calculation that is in 
parallel with the results in Fig. [3] but takes into account 
static disorder with the noise amplitude A = 0.5. As 
demonstrated in Fig. [71 in the presence of such a high 
noise level, the quantum population of spin excitation is 
still successfully transferred to the last spin of the chain 
(peak probability around 99%), with the excitation pro- 
file almost unaltered as compared with the noiseless case 
studied in Fig. [3J For an even higher fluctuation level, 
e.g., A = 0.7, the spin excitation profile is seen to grad- 
ually disperse as it is transported along the chain. 

Figure [H] displays results for an initial Gaussian excita- 
tion profile also considered in Fig. [4J but in the presence 
of static disorder characterized by A = 0.5. We find 
that the severe disorder can slightly change the shape of 
the spin excitation profile (though not so evident in Fig. 
8) and hence the peak value of the probability profile 
slightly fluctuates during the controlled transport. Note 
however, the area enclosed by the main probability pro- 
file is found to be around 0.99 at all times. At time 
t = 21500, the peak probability of the spin excitation 



FIG. 8: Adiabatic transfer of an initial Gaussian profile of 
spin excitation along a chain of 101 spins, in the presence of 
static disorder characterized by A = 0.5. Other parameters 
are the same as in Fig. [4] Panels (a), (b) and (c) are for 
t = 6000, 12000 and 21500. Results here are very similar to 
those shown in Fig. [4] 



profile, still as large as 0.99, has been transferred to the 
100th site as in the noiseless case of Fig. [4] In another 
sampling case for aA = 0.7, the peak excitation prob- 
ability that is transferred to the last spin decreases to 
0.94. All these results clearly demonstrate the robust- 
ness of our adiabatic transport scheme to high-level static 
disorder. 

We have also examined the robustness of our adia- 
batic transport scheme to dynamic disorder. To model 
time-dependent fluctuations in the spin-spin coupling 
strength, we now let each S n be given by the sum of 
ten oscillating functions, i.e., S n = J2i=i Acos(ujit + tfii), 
where u>i are random frequencies distributed in [0,w max ], 
and 4>i are random phases uniformly distributed in [0, 2n] . 

Interestingly, our numerical experiments indicate that 
effects of dynamical disorder modeled above depend 
strongly on cj max , i.e., the cut-off frequency of the 
dynamic fluctuations. Introducing disorder more fre- 
quently, i.e., introducing a larger w max , can lead to much 
decreased peak probability transferred to the last spin. 
In particular, we find that for w max < 0.1, the effects of 
the dynamic disorder are essentially analogous to what 
is found for static disorder. For larger w max , the dete- 
rioration of the adiabatic population transfer becomes 
considerable for the same noise amplitude A. Figure [5] 
displays the results for A = 0.025 and ui max = 0.1. Note 
that for A = 0.025, the amplitude of the noise is very 
large because the total fluctuation is a sum of ten func- 
tions oscillating at the same amplitude A. It is seen from 
Fig. [9] that for such a case of dynamic disorder, the spin 
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FIG. 9: Adiabatic transport of spin excitation for an initial 
state localized exclusively at n — 0th site along a chain of 101 
spins, in the presence of dynamic disorder characterized by 
A — 0.025 and U) maK = 0.1. Other parameters are the same 
as in Fig. [7] Panels (a), (b), and (c) are for t = 7000, 14000 
and 20000. 



excitation travels almost unaffected along the chain, thus 
confirming again the robustness of our adiabatic scheme. 
However, upon an increase in w max , e.g., u; max = 1.0 (so 
noise frequency becomes comparable to the characteristic 
coupling strength J), the dispersion of the spin excita- 
tion profile becomes evident in Fig. [TO] The situation 
can be certainly much improved if the noise amplitude A 
is decreased. Because similar results are also found for 
Gaussian excitation profile as initial states, we conclude 
that the noise spectrum of dynamic disorder can play an 
important role in affecting the robustness of our adiabatic 
transport scheme, especially when the noise amplitude is 
very large. 

Although our adiabatic scheme is seen to be robust, 
it can be expected that the very existence of disorder 
should limit the threshold speed of the moving potential. 
Put alternatively, for a larger moving speed (which sat- 
isfies the adiabatic condition less), the robustness of our 
control scheme to disorder is expected to decrease. This 
trend is indeed found in our numerical experiments. To 
characterize precisely how an increasing moving speed of 
the control potential affects the robustness is certainly 
beyond the scope of this work. 



V. ADIABATIC TRANSPORT IN A DUAL SPIN 
CHAIN 

As demonstrated in previous sections, our adiabatic 
scheme based upon a moving parabolic potential offers 



FIG. 10: Adiabatic transfer of spin excitation for an intial 
state exclusively localized at the n = 0th site along a chain of 
101 spins. aj max = 1.0, and other parameters are the same as 
in Fig. HD 



a simple and robust approach to transferring quantum 
population along a spin chain. This scheme requires a 
strong parabolic field if it is globally parabolic, but even 
this requirement can be greatly weakened if the initial 
spin excitation profile spans a few sites. Further, one 
does not really need a globally parabolic field to realize 
this adiabatic scheme: it suffices for the parabolic field 
profile to be wider than the spin excitation profile. With 
these considerations we may argue that the well-known 
dispersion issue in quantum information transfer along 
spin chains is essentially solved by our adiabatic scheme. 
Nevertheless, as also mentioned earlier, one important is- 
sue still remains open. That is, for the sake of quantum 
information transfer, how to take care of the quantum 
phase of a quantum state to be transported? Indeed, a 
moving external potential induces extra dynamical phase 
to the spin chain, and such a dynamical phase depends 
on the details of the control potential. These facts mo- 
tivate us to seek an encoding approach that can protect 
a quantum state from the additional dynamical phases 
induced by the moving potential. 

Fortunately, the idea of using a dual spin chain, first 
proposed to overcome the disorder and dispersion issues 
in quantum information transfer in spin chains [i~7l [2fj| , 
offers a promising solution. Specifically, we propose to 
combine our adiabatic transport scheme with the dual 
spin chain scheme. Then, because each individual sub- 
chain acquires identical dynamical phases from the same 
external moving potential, the relative quantum phase 
between the two sub-chains is certain, and hence quan- 
tum information encoded in the dual spin chain can be 
transported without suffering from the extra uncertain 
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dynamical phases. 

Consider then a quantum channel consisting of two 
identical parallel spin chains subject to the same external 
parabolic potential, 



J 



N-l 



n=0 



A' 



J2 C 



(n - n ) 2 W| 



(9) 



where i = 1,2 indices label the two sub-chains. Suppose 
the quantum state to be transferred is given by |4>) = 
a|0)+/3|l). Such a state can be encoded into the quantum 
channel prepared in the following entangled state, 



l*(0)> - a\ S ) 



(i) 
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as a superposition of two components: the n — 0th spin 
in the second (first) sub-chain being flipped and the first 
(second) sub-chain in its ground state denoted by |g). 
Note that for each sub-chain at most one spin is flipped 
and the associated dynamics will be restricted to the 
ground state or the subspace of one flipped spin. 

This encoding can be extended to cases of entan- 
gled Gaussian wavepackets in a straightforward manner. 
However, for convenience here we discuss only cases aris- 
ing from the initial state given by Eq. (jTUJ) . After the 
independent evolution of the two sub-chains for a total 
duration of r under the action of the moving parabolic 
potential, the quantum state of the dual spin chain is 
given by 



l*(r)> = $>(t)|*»> 



(11) 



where = a|g) (1) <g)|n)( 2 )+/3|n)W<g)|g)( 2 ). Evidently, 
though each c„ contains the extra quantum phases in- 
duced by the external moving potential, this factor is 
identical for the two state components of |3? n ). As al- 
ready demonstrated in our numerical experiments using 
a single spin chain, the profile of |c„| 2 should also be 
highly localized, and the time of arrival of the peak value 



of 



at the last spin can also be directly calculated 



from the moving speed of the parabolic potential. 

Analogous to the original dual spin chain scheme, at 
the end of the adiabatic quantum transport the final state 
|\I/(t)) can be decoded by applying a CNOT operation 
to the last two Nth spins of the dual chain. Upon this 
operation the final state is transformed to 

N-l 

c„(r)|$„) + c n (t) [a|g)« + /3|N)«] ® |N)( 2 >(12) 

n=Q 

As such, by measuring the last spin of the second sub- 
chain, one gains important information about the trans- 
port. In particular, if the measurement outcome is spin 



up, then the initial state |$) = a\0) + /3\1) has been suc- 
cessfully transferred to the last spin of the first sub-chain, 
with probability |cat(t)| 2 ; if the outcome is spin down, 
then the quantum state transfer is unsuccessful and one 
needs to wait for more time to perform additional mea- 
surements. 

Significantly, because our adiabatic population transfer 
scheme can ensure a very high probability of excitation 
transfer to the last spin, the probability of spin-up mea- 
surements can be guaranteed to be very high (arbitrarily 
high if there were no restriction on the field strength). 
This hence overcomes, at least theoretically, one main 
disadvantage of previous dual spin chain schemes where 
too many measurements may be required for high fidelity 
quantum state transfer. Further, at the end of the adi- 
abatic transport, the spin excitation is automatically lo- 
calized at very few end spins. So it also becomes un- 
necessary to perform fast measurements at a particular 
time. Instead, one can choose measurement times at will 
so long as the moving parabolic potential has reached the 
last site of the spin chain. This makes it clear that our 
adiabatic scheme, when combined with quantum phase 
encoding schemes, can find important applications in 
quantum information transfer (in addition to quantum 
population transfer). 



VI. CONCLUSIONS 

In this work we have presented a simple and robust 
scheme to realize adiabatic population transfer in spin 
chains. The additional resource needed is a slowly mov- 
ing external parabolic magnetic field. The basic mecha- 
nism is the adiabatic following of a quantum state with 
the movement of a separatrix structure in the classical 
phase space of a pendulum analogy. In particular, we 
have shown that our scheme can be used to transfer spin 
excitation from one end of a spin chain to the other end, 
with the initial excitation profile being a localized trun- 
cated Gaussian wavepacket or exclusively localized at a 
single spin site. It is found that much weaker external 
field is needed for adiabatic population transfer if the 
initial excitation profile covers a few spin sites. Effects 
of static and dynamical fluctuations in the spin-spin cou- 
pling strength are also computationally studied, confirm- 
ing the robustness of our adiabatic population transfer 
scheme. Realizing the robust population transfer with 
small dispersion, we have also proposed to apply our ap- 
proach to a dual spin chain such that robust quantum 
information transfer can be realized with important ad- 
vantages. We hope that our theoretical scheme can moti- 
vate experiments using various implementations of a spin 
chain Hamiltonian, such as cold atoms in an optical lat- 
tice and electron tunneling in an array of quantum dots. 

The central idea of this work, namely, using a slowly 
moving external potential to adiabatically transfer spin 
excitation, might be useful for other applications as well. 
For example, one may consider distributing entanglement 
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along a long spin chain in a controlled fashion, by use of a 
control potential that has two components slowly moving 
in opposite directions. Another interesting application is 
related to studies of quantum signal amplification with 
spin chain models. Recently, an interesting connection 
between quantum state transfer and quantum state am- 
plification is revealed [27j |. In this regard, our adiabatic 
scheme might also help design a new and useful approach 
to controlled quantum amplification using a slowly mov- 
ing external potential. 
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